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. Abstract 

^\ . ^_ . . 

The integrated branching fraction of the process B — > Xs I I is dominated by 

resonance background from narrow charmonium states, such as B ^ Xgip — > 
Xg , which exceeds the non-resonant charm-loop contribution by two orders of 
magnitude. The origin of this fact is discussed in view of the general expectation 
of quark-hadron duality. The situation in i? ^ Xg l^l~ is contrasted with charm- 
penguin amplitudes in two-body hadronic B decays of the type B — > tttt, for which 
it is demonstrated that resonance effects and the potentially non-perturbative cc 
threshold region do not invalidate the standard picture of QCD factorization. This 
holds irrespective of whether the charm quark is treated as a light or a heavy quark. 



1 Introduction 



It is a well-known fact that the resonant transition B — > Xgip — >■ Xs I'^l" exceeds the 
non-resonant short-distance process B ^ Xs by about two orders of magnitude, 
even when considering the rate fully integrated over the dilepton invariant mass. Since 
the process is inclusive over the hadronic final states, this is sometimes interpreted as 
a particularly striking failure of global parton-hadron duality, which says that the sum 
over the hadronic final states, including resonances, should be well approximated by a 
quark- level calculation [1]. The origin of the problem is apparently rooted in the fact 
that the process B ^ Xg I'^l" contains a penguin contribution with a charm-quark loop, 
which has only a small effect on the partonic calculation but leads to a large resonant 
enhancement in the real world. Charm-quark loops are relevant in other contexts as well, 
such as in the calculation of the charm production cross section in e^e~ annihilation from 
the imaginary part of a current correlator, and in the calculation of penguin amplitudes 
in non-leptonic B decays in the QCD factorization framework [2, 3, 4]. Should one 
expect large violations of global parton-hadron duality due to charmonium resonances 
in these cases as well? The purpose of this note is to explain the origin of the resonance 
dominance 'm. B ^ Xs I'^l' decays. We then investigate the differences between this 
process, where duality should not be expected to hold, and non-leptonic B decays, where 
we demonstrate that it holds in the heavy-quark limit. 

The contributions of charmonium states to the different observables studied in this 
note are described in terms of resonance contributions to a correlation function n(g^). 
The question of when large deviations from quark-hadron duality arise is tightly linked 
to how this function enters the formulae for the various decay rates. After a brief re- 
view of some basic facts about B ^ Xg decays, we study the correlator of two 
charm-quark currents using the operator product expansion (OPE) as well as a hadronic 
picture in terms of Coulomb resonance states. We then construct a toy model, which 
elucidates important features of the correlator and their relation to violations of quark- 
hadron duality. Returning to the case of -B — > decays, we derive numerical 
estimates explaining the observed large violation of duality in a semi-quantitative way. 
Finally, we consider the question of duality violations in charm-loop penguin contribu- 
tions to exclusive hadronic decays such a.s B ^ titt. Our discussion in this note focuses 
on the inclusive leptonic process B ^ Xg , for which duality violations in the charm 
resonance region are particularly pronounced and experimentally well studied. We em- 
phasize that an analogous discussion holds for the related radiative process B — > Xg'-f, 
for which the resonant transition B Xgip followed hj ^ X7 also exceeds the total 
inclusive short-distance contribution. 

2 Brief review of B — > Xg decays 

The total inclusive branching fraction for the process B ^ Xg l^l' may be written as [5] 
B{B Xg /+/-)sD = rs ^ ^ \V,g\' + \Cio? ^ . , 
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C9,Cio 12 
Figure 1: Short- distance (left) and charm-penguin contributions (right) to 6 — > sl'^l". 



where we have neglected the contribution from the magnetic dipole operator Qj^, which 
is immaterial to our discussion. The subscript "SD" indicates that this expression is the 
short- distance contribution to the branching fraction, computed using partonic matrix 
elements of the effective weak Hamiltonian. The quantity (jCgp) is understood to be the 
effective coefficient |Cg^p, which includes the penguin-type matrix elements, averaged 
over the dilepton-mass squared: 

(ICgH = 2 f'ds (1 - sY (1 + 2s) Cg^D^ + 3a2 h{z, s) ' , (2) 
Jo 

where h{z,s) with z = rric/mh and s = q'^/ml is the charm-penguin loop function, 
and 02 = C2 + Ci/3 ~ 0.12 [5].^ The leading contributions to the decay rate are a 
large logarithmic term ~ ln(Mi^/mb) in Cg from the ultraviolet contribution to the 
charm-penguin amplitude (high-energy scales between rrib and Mw) and the top-quark 
contributions to Cg and Ciq. Since C^°^ ^ —Ciq ~ 4 and h{z,s) is of order 1, we 
conclude that the low-energy contribution (scales at or below nih) of the charm penguin is 
subleading. Numerically, ((|Cgp) + |Cion/2 f» 16 in (1). We also use a = a{mh) = 1/132 
and rrih = 4.8 GeV. Feynman diagrams illustrating the short-distance and charm-loop 
contributions are shown in Figure 1. 

On the other hand, the branching fraction for B Xgip reads (with r = M^/ml) 

B{B Xs^) = TB \VM' 4 (1 - (1 + 2r) ^ 1.7 ■ 10"^ (3) 

This falls short of the experimental value B{B Xgip) = (7.8 ± 0.4) ■ 10"'^ [6], since 
the inclusive ip production rate is probably dominated by colour-octet production [7, 8]. 
Phenomenologically, the experimental rate may be accounted for by introducing a factor 
K = 2.16 and replacing 

a2 K,a2 ■ (4) 

Multiplying the experimental value with the leptonic branching fraction B(?/^ — ^ = 
0.06, we obtain the resonance contribution B{B — * Xsip /"*"/") ^ 4.7 ■ 10"^ and 

thus 

^V- = p/D . V ~ ^0 . (5) 



B{B ^Xsl+l- 



jSD 



^Our cocfRcients Ci 2 are called C2 1 in this reference. 



2 



This estimate is only slightly modified by contributions from the higher charmonium 
resonances. Since the denominator includes the partonic charm-loop contribution, we 
conclude that the quark-level calculation of the total branching fraction yields a result 
that is far smaller than the contributions from the charmonium resonances alone, indi- 
cating a gross failure of global quark-hadron duality. 

In order to include the resonance contributions in a phenomenological description of 
the global features of the B ^ Xg spectrum, Kriiger and Sehgal (KS) have suggested 
to relate the charm-penguin loop to the hadronic cc vacuum polarization function 11 (g^) 
using a dispersion relation [9]. This assumes that QCD interactions between the charm 
loop and the remaining quarks in the process can be neglected. On the other hand, 
the KS model offers the possibility to take resonances and other hadronic contributions 
related to the charm loop into account. Employing the dispersion relation gives the 
correct scheme dependence of the loop function and avoids a double counting of partonic 
and hadronic degrees of freedom. The ansatz of [9] amounts to the replacement of the 
penguin function h{z, s) by 



The vacuum polarization function n(g^) will be defined in (7) below. Its imaginary part 
can be obtained from experimental data and related to n(g^) through the dispersion 
relation (8). An explicit parametrization with six Breit-Wigner resonances and a con- 
tinuum contribution can be found in [9]. To account for the experimental B — > X^ip 
rate discussed after (3), KS include the phenomenological enhancement factor k in (4), 
which increases the ip resonance term. The global spectrum for B ^ Xs decays 
with resonances included using the KS approach is shown in Figure 2. 

From the preceding discussion it is clear that important information on the hadronic 
dynamics of B ^ Xs is contained in the vacuum polarization function, which we 
briefiy review in the following section. 

3 Vector-current correlator 
3.1 Parton-hadron duality 

In discussing the validity of a perturbative quark-level computation of the penguin am- 
plitude, the concept of quark-hadron duality is crucial. We will therefore recall basic 
aspects of duality for the correlator of two vector currents. This is a particularly simple 
situation, but it is also relevant for the penguin amplitude in the approximation where 
further gluon attachments to the charm loop are neglected. 
The correlator is defined by 




4 IQtt^ r 



n(g2) _ n(o) 



(6) 
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(7) 



3 



2.5 




Figure 2: Differential B — > branching fraction as a function of s = q'^/ml = 

mf+i-/ml, including the effect of charm resonances in the KS method (solid line). A 
phenomenological factor of k = 2.16 is used for the resonance amplitude. For com- 
parison, the dashed curve shows the same quantity obtained within a purely partonic 
calculation. The short-distance contributions are computed at next-to-leading logarith- 
mic order. 



where = C7^c for the charm quark, and c ^ u,d,s for the case of light quarks. The 
scalar function 11 (g^) obeys the dispersion relation 



n(g^)-n(0) = ^ T^7=^, (8) 



TT Jo t t — — ie 

which is defined in the Euclidean domain < and may be analytically continued to 
the time-like region > 0. The expression on the left-hand side of (8) can be computed 
in QCD by means of an OPE as a power series in l/g^, with the coefficient of each term 
given as an expansion in the strong coupling a^. We shall denote this quark-level (OPE) 
result for n(g^) as Uglq"^). To lowest order in the OPE, and omitting 0{as) corrections, 
the correlator is given by 

Iiq{q^) - ng(0) = -— / duu{l - u) In ^ ^ -TTTT ' 



(9) 

with m the quark mass and = 3 the number of colours. The second expression in (9) 
is the asymptotic behaviour for large g^. The imaginary part is non-zero only if q^ > Am? 
and reads 

imn,(g2) = -^\i--r[^ + -r]- (lo) 
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On the other hand, eq. (8) is vahd for the exact, hadronic correlator Uhiq^). The 
imaginary part of Uh is related to the observable R ratio through^ 



Rc 



aye^e 



hadrons (cc)) 



(11) 



(T(e+e~ 



where Cc = 2/3 is the charm-quark charge. Using (8), Ilhiff) may then, in principle, be 
obtained from experimental data. 

To any finite order in ag and the OPE quantity I^q{q^) is expected to approx- 
imate the true function 'n.h{q^) in the sense of an asymptotic expansion. For < 0, 
with sufficiently large the difference between Uq and Uh should then be systemati- 
cally reduced with the inclusion of more terms in the OPE, up to a residual uncertainty, 
which is expected to be exponentially suppressed as exp(— c Ig^l/Aqcj-,) with an 0{1) 
constant c. For > the deviation of the smooth function Uglq"^) from n/j(g^) might 
also involve oscillatory terms in g^. The approximation of Uh by Ilg will then in general 
not work point by point, but only if > is large enough or Imll is averaged over 
a sufficiently wide range in q^ [1]. The correspondence between Uq and point by 
point in is referred to as local quark-hadron duality, the correspondence of averages as 
global quark-hadron duality. An illuminating discussion of these concepts can be found 
in [10, 11]. 

3.2 Vector-current correlator in the Coulombic limit 

The concept of duality is of particular importance for light-quark correlators, where the 
resonance region is always non-perturbative. For massive quarks, we may also consider 
the opposite limit, where not only 2m, but also the binding energy of the quarkonium 
resonances is much larger than the strong-interaction scale. This limit is referred to 
as "Coulombic", since the dominant binding force is then the colour-Coulomb force. 
Although charmonium is not a Coulomb bound state, it is instructive to discuss the 
Coulombic limit of the charm-quark vector-current correlator in the present context. 
Some of the results of this section will also be used in the numerical estimates below. 

Near ^ 4m^ the spectral function lmllh{q^) receives prominent contributions from 
a series of cc resonances (and other hadronic states with open charm), which cannot be 
captured by a fixed-order perturbative expansion. The exchange of Coulomb gluons gives 
rise to corrections of order Og/v, where v is the velocity of the heavy quarks [12]. In the 
threshold region i; ~ the entire ladder of Coulomb gluons needs to be resummed in 
order to account for the detailed behaviour of n(g^). When this is done, one finds for 
the correlator after MS subtractions [13] 



^The relation is approximate, since cc pairs may annihilate or be produced from gluon splitting. 



asCpN 1 
871 X 




1 , -Am^E 1 



- + 7£; + ^(1 - A) 



(12) 
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with 

A = ^y^, E = ^^-2m^ + ^e. (13) 

Here il){z) = d\nr{z)/dz denotes the digamma function, and Cp = (iV^ — 1)/{2N). 
Expression (12), vahd for ^ 4m^, contains poles for positive integer values X = n, 
which correspond to the charmonium binding energies 

in the colour- Coulomb potential. The residues of these poles determine the decay con- 
stants /„ of the S-wave bound states at energy En. The pole part of (12) is thus equivalent 
to 

oo r2 

where M„ = 2mc + En and = f^/n^- The decay constant of the ground state is 
given in terms of the wave function at the origin by 

,2 _ 4JV|V(0)P _ 16 2 3 |,,„M2_ 1 /'"c 



fi-^i^-9i<'''' \m\'--[^C,a.) . (16) 

The validity of local vs. global parton-hadron duality in the Coulombic limit is related 
to the issue of perturbative resummation. When sufficiently wide and smooth averages of 
Im n or n over the resonance region are taken, global duality holds in the sense that these 
quantities are well approximated by fixed-order perturbative and OPE computations 
without any need to sum the Coulomb ladders. In addition, in the Coulombic limit one 
may expect that local duality holds as well, since the physics of the threshold region 
is perturbative. It is evident that a local correspondence of Imll/i and Imll^ requires 
the perturbative resummation of Coulomb gluons discussed above. However, even in 
the ultra-heavy-quark limit there is a non-perturbative contribution from the strong- 
interaction scale to the charmonium binding energy [14] . Thus, for local duality to hold 
this contribution must be small compared to the width of the resonance. 

The widths of the resonances are not present in (12) and (15), since in the Coulombic 
limit they are of higher order in small coupling constants. A further resummation is 
required that replaces — > — iMnTn in the denominator of (15). Concentrating 
on n = 1, the resonance can decay via the strong interaction into a three-gluon final 
state or into a fermion pair (e, /i, u, d, s) through a virtual photon. The total decay rate 
of the ground state is then given by 

= r(^ ^ 3^7) + 4r(^ ^ e+e") , (17) 

where 



r(^^3«) = 12(^!^A, r(.^^e-e-) = i^. (18) 

^' 243 ' ^ 27M^ ^ ' 
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This expression requires only that the charm-quark mass is large compared to the strong- 
interaction scale. In the Coulombic limit, one may further use the result (16) for 
Note that the second term in (17) constitutes a sizable fraction of about 20% of the total 
rate, despite the electromagnetic origin of this contribution. 



4 A toy model 



Before we return to a discussion oi B ^ Xs decays, we shall consider a toy model 
in which the hadronic part of the amplitude is exactly the current correlator discussed 
above, such that the role of resonances and quark-hadron duality is exhibited in a par- 
ticularly transparent way. To this end, let us hypothetically assume the existence of two 
"leptons", /i with a large mass mi and I2 with mass 1112 = 0, and the effective weak 
Hamiltonian 



n 



G r.. 



V2 



{hh)v-A icc)v-A - {khW-A itt)v-A ■ (19) 



Since we are interested in the QCD dynamics of the decay, the flavour aspects of the 
model are unimportant for our discussion. All particles are assumed to have standard 
strong and electromagnetic interactions. Then Ti-eff gives rise to a loop- induced process 
li —>■ I2 e^e~ via charm- and top-quark penguin diagrams with a GIM-like cancellation 
between them. The corresponding decay amplitude reads 

A{h ^ I2 e+e-) = e,e2 Ii{q^) l2Y{l - l,)h e^^e . (20) 

Here 11 = lie — lit is given as the difference between the charm and top contributions. 
We take rrit > mi, and thus Imll comes only from the charm sector. The correlator 
n(g^) fulfills the dispersion relation (8), where 11(0) is fixed by our model and can be 
computed in perturbation theory. To leading order one finds 

TV m^ 

n(o)^n,(o)-n,(o) = -— in^. (21) 

livr^ m^ 

The form of (20) for the amplitude holds to lowest order in G and e^, but to all orders 
in the strong coupling. In fact, since the quark loops have no QCD interactions with the 
other, purely leptonic parts of the amplitude, the quantity 11 in (20) can be considered 
as the exact hadronic correlator in QCD. From (20) we obtain the differential decay rate 
(with s = q^/ml) 

dnh-^he^e~) ^ GW, _ ^ 2s) \Ili<f)\\ (22) 
ds 2 '~l IX 

We shall now investigate to what extent the hadronic function n(g^) in (22) may 
be approximated by a quark-level calculation. This is clearly the case for values of 
at which local duality is a reasonable approximation, that is if q^ ^ 4m^ is well above 
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the narrow-resonance region, or if < 4m^ lies sufficiently far below the charmonium 
threshold (we will make the latter statement more precise at the end of this section). 
However, in a region of where prominent resonances exist, local duality is certainly 
badly violated. Then |n(g^)p and also Imn(g^) cannot be computed locally in the OPE 
approach. For Imn(g^) the agreement of the full hadronic correlator with the quark- 
level calculation can be improved by integrating with a smooth weight function over a 
sufficiently large interval in (global duality). However, in (22) the correlator is squared 
before the g^-integral can be performed to obtain the total rate for li — >• l2e'^e~. As a 
consequence, the impact of sharp resonances is not averaged out even for the integrated 
rate. This is in fact our main point: The concept of global parton-hadron duality should 
not be expected to apply to averages of the quantity |H(g^)p. 

To make this point clearer, we consider the contribution of a single ip resonance to 
H(g^), omitting for the moment other resonances and non-resonant terms. The correlator 
then has the approximate form 



q^-Mj+iM^T^ 



n(g^) = - „. .... , (23) 



where M^,, and denote the decay constant, mass, and total width of the resonance.^ 
It follows that 

imntf) = -^4J#IW. (24) 



- AQ)2 + Mm 



while 



|n(g2)P = — ^§ ^ = _S_iniH(g2). (25) 

It is important to note that in the small-width limit |H(g^)p is more singular than 
ImH(g2). 

Concentrating ffist on the absorptive part of the correlator, we note that it can 
deviate by a large factor from the partonic result. For instance, on resonance we have 



ImH(M^) = ^ ^ ImH,(M:^) = 0{l) . (26) 



r2\ _ 



At the V resonance, for which = 3.097 GeV, = 0.401 GeV, and F^ = 93.2 ■ 
10-6 GeV, 



M,,T 



560 , (27) 



■^If n((7^) is calculated within QCD. then strictly speaking excludes the electromagnetic contribu- 
tion in (17) to the total decay width. In the following, it will be useful to consider as an adjustable 
parameter. For instance, we might consider a variant of the model, in which Tl{q'^) is the correlator of 
a current with different quark flavours, so that the QCD contribution to the width of the resonance 
is exactly zero. In this case F^ is proportional to the arbitrary coupling constant of a hypothetical 
interaction by which the resonance decays to e+e^. 
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corresponding to a huge violation of local duality. This behaviour has both a parametric 
and a numerical origin, as can be seen from (18), which gives, using as{M^) = 0.25, 



•^^ ~ ZZr ~ A^o (2R) 

M^T^ 40(7r2-9)«3(M^) ' ^ ) 

in reasonable agreement with the experimental result (27). In spite of this, the resonance 
contribution to the integral 
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dqHmUiq^) ^7ifl<^ml (29) 



is parametrically small compared to the integral over the non-resonant contribution, 
which is of order ml when ml is large. This is consistent with global duality. In fact, 
the integral over the partonic expression for Imn(g^) is dual to the integral over the 
hadronic expression including the contribution from the resonances. This property is 
used in the determination of the charm-quark mass from QCD sum rules [15, 16], as well 
as in the determination of from r decays [17], in which case it is applied to light-quark 
resonances. 

We now contrast this to the situation for the square of the correlator. For the narrow 
ijj resonance, using (25) and 

lmU{q')^7rf}6{q'-Mi), (30) 
we obtain for the resonance contribution to the integral 

rvw)r-7r4x-^. (31) 

This should be compared with the non-resonant contribution, which remains of order 
mf. The appearance of the factor /^/(M^r^) with the width in the denominator allows 
the resonance contribution to be arbitrarily large. Since expression (20) is exact in QCD, 
it is clear that global duality cannot possibly hold in this case and may be violated by 
an arbitrarily large amount (to the extent that we can adjust the resonance width at 
will). Substituting the resonance expression for |n(g^)p into (22), we find 

rrcs('i '2 e^e )= / ds = T{li ^ hip) , (32) 



ds T. 



where (with r = Ml/m\) 



r(/i - h^) = (1 - rf (1 + 2r) . (33) 

This reproduces the result for the rate of the decay chain li —>■ I21I) — I2 e^e~ in the 
narrow-width approximation. 
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The ratio of the resonance rate to the non-resonant partonic rate can now be 
estimated as follows. First, the partonic rate is obtained from (22) using the partonic 
expression for n(g^). In the one-loop approximation we find 

IGtt^^ / 2n 8, mt , 20 , 4 2 , , ; f 2arctan^^; x > 1 



.-n^fe^'-'-'^^^-^^T-^r-if^-VU-li,^^,^^^ ^^^^ (34) 



where x = Aml/q'^. For typical values of the parameters (e.g. = 1.4 GeV, rrii = 

rrib = 4.8 GeV, rrit = 167 GeV) the minimum of |nq(g^)| is located at = 0, while its 
maximum occurs at = 4m^. At these values of 

= 7^ (^^^ + c] , (35) 



2n'^ \ rric 

with c = at = and c = 4/3 at g^ = 4m^. The total rate is well represented 
by treating |ng(g^)| as constant with an intermediate value of c = 1/2. It can then be 
estimated to be 

G'^a'^ml f, nit 1~ 



r(,.^,,e-e-)3o«^(ln^ + -) . (36) 



The resonant contribution is 



r.es(/i - - k e+e") = ^y,ff^' (1 - r)2 (1 + 2r) . (37) 



27M/,r 



We thus find for their ratio 



R _ r.e.(/i^/2^-^/2e+e-) _ 87r^(l-r)^(l + 2r) 4 /| 

T{h - h e+e-)sD (ln(m,/m,) + 1/2)^ m? M^F^ ^ ' ^''''^ 

In this example the first factor is about 55, f^/m1 = 0.007, and the last factor has 
been given in (27). We observe that the moderate suppression of the channel /i — >■ ^ip 
by 55/|/mf ~ 0.38 is overcompensated by the enhancement factor /|/(M^F^), whose 
origin we have discussed above in the context of (25). Since the latter factor is very large 
because of the narrow width of the ip resonance, we have Fj-cs ^ Fs^. As anticipated, 
this result excludes the applicability of the quark-hadron duality hypothesis. 

The deeper origin for the absence of global duality can be seen in the fact that there is 
no OPE for the total rate of li — * I2 e+e" decay. Constructing the rate as the absorptive 
part of the li li forward-scattering amplitude, 

{h\lmi I d^xHeffix) Heffmh) , (39) 

we would need contributions with an e~^e~ loop, coupled to charm- and top-quarks by 
photons, and we would have to require that the cut passes through this loop in order to 
obtain the rate for li I2 e~^e~ (see the first graph in Figure 3). This does not correspond 
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Figure 3: Cuts through the h li forward-scattering diagrams that contribute (a) to 
the li — s> I2 e^e~ decay rate and (b) to the inclusive hadronic decay rate for li ^ I2 + X. 
See text for further explanation. 



to an OPE of (39), which would imply a sum over all cuts, not just a restricted set. In 
this way r(/i — > l2e^e~) is seen not to be a truly inclusive quantity, for which global 
duality would be expected to hold. Rather, the selection of a final state with an e^e~ 
pair represents a more "exclusive" choice, even for the integrated li I2 e^e~ rate, 
which leads to an integral over |n(g^)p. The situation encountered here is similar to 
that for the radiative decay B — »• ^^7. As emphasized in [18], in this case contributions 
to the decay rate for which the photon is not part of the local operators in the effective 
weak Hamiltonian cannot be obtained from an OPE. They correspond to a subset of 
cuts analogous to those in Figure 3. 

On the other hand, a very different situation occurs for the inclusive hadronic decay 
/i I2X (see the second graph in Figure 3). In that case no restriction is placed 
on the cuts, and an OPE can be apphed to (39). Similarly to the rate of hadronic r 
decay [17], the decay rate is given by a weighted integral over the imaginary parts of the 
vector- current (11) and axial- vector current correlators (11^'^), 

T{h - /2X) = — ^ / dsil- sf \{l + 2s) (lmn(g2) + lmll\{q')) + Imni(g2)" , 

(40) 

and global duality works in the same way as for the charm contribution to the e'^e — >■ 
hadrons cross section. 

Finally, we briefly return to the low-g^ region in li I2 e^e~ decays mentioned after 
(22). Here a quark-level calculation is justified if is sufficiently below the if) resonance. 
We estimate Q'^ax' maximum value of g^, up to which a quark- level calculation of 
n(g^) can be trusted. For close to zero the quark picture is reliable, and the if)- 
resonance contribution is only a small part in a hadronic representation of 11. As q^ 
gets close to M|, the resonance contribution dominates the correlator while the partonic 
result for 11 is too small. Therefore, as an estimate for q^^-^ we use the point in q^ where 
the one- loop partonic result equals the ■^/'-resonance contribution to 11 (g^) — 11(0). Using 
[n(g^) — n(0)]^ = (/|/M|) g^/(M| — q^) and the partonic expression in (34) close to 
threshold q^ = Ami ~ "we obtain g^^^^ = M| — 37r^/|/2 ^ 7GeV^. Here we have 
considered the limit /|/M| <^ 1 in order to obtain a simple analytic expression. We then 
have H(g^jj^^) — H(0) = 0.04 for the one-loop result, 0.05 for the ■i/'-resonance contribution. 
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and 0.06 for the full hadronic expression in the KS representation. Since n(g^) — 11(0) is 
subdominant in comparison to 11(0), given in (21), the quark-level estimate of Il{q'^^^) 
is still rather accurate. A similar discussion applies to B ^ Xg . 



5 Charm resonances 'u\ B — > Xg l^l decays 

We now return to the discussion of -B — Xgl^l^, following the analysis of the toy 
model considered in the previous section. The branching fraction for the decay chain 
B Xgip Xg is obtained by multiplying (3) with 



where the total rate is given in (17). Replacing 02 — ~ 0.26 in (3) and using 
(1), we then have for the ratio of the resonant and the partonic rate 

_ B{B ^ Xsj, ^ Xg in-) _ 5127r5/,^ (1 - rf (1 + 2r) /|_ /| 

- B(i? ^X,/+/-)sD 9((|C9|2) + ICloP) ''ml'' M^V^ ' ^ ' 

The first factor is about 23 numerically. The first two factors give approximately 0.16. 
The large enhancement from /^/(M^r^) = 560 overcomes the suppression ~ 0.16, and 
we recover pa 90. This explains the size of R^ already quoted in (5). 

It is interesting to consider the heavy-quark limit mfe, rric ^ Aqcd with mc/mf, fixed, 
where the ip resonance is asymptotically a Coulombic bound state. If we still assume 
a <^ as, then (16) and (18) imply 

We see that formally R^ = 0{1) in the heavy-quark limit. However, the expression (43) 
contains a large enhancement from numerical factors and from the running of ag between 
the momentum scale of the charmonium bound state, entering and the scale M^, 
which is relevant for F^. Using as{mcv) / as{M^) ^ 2, the asymptotic formula (43) gives 
R^ ~ 60 and thus reproduces the bulk of the resonance enhancement. In the heavy- 
quark limit a resummation of Coulomb ladders would automatically account for the large 
■j/'-resonance contribution in the total rate oi B ^ Xsl^l~ . Upon Dyson resummation 
of finite-width effects into the resonance propagator, expression (43) is obtained. 

The large charm-loop effect may be contrasted with the case of a light resonance p, 
where Fp, Mp, /p ~ Aqcd, and therefore the ratio i?p, defined in analogy to (42), is 
strongly suppressed. Estimating the first factor to be about 10-50 and using /^/(MpFp) 
= 0.38, one finds 

Rp ^ [10-50] X ^ X ^ [0.007-0.036] . (44) 
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Figure 4: Penguin diagram with a charm-quark loop contributing to B — » M1M2 decays. 
The curly line may either represent a gluon in the case of a QCD penguin graph, or a 
photon for an electromagnetic penguin. 

Even though we cannot expect quark-hadron duality to hold, in practice the p-resonance 
contribution is negligibly small. This scenario is relevant for up-quark penguins in S — ^ 
Xs,dl^l~ decays. 

6 Penguins with charm in B — > tttt decays 

It is interesting to contrast the situation in the rare leptonic process B ^ Xs with 
the penguin amplitudes in exclusive hadronic decays such as S ^ tttt, where similar 
diagrams with intermediate charm quarks contribute [19, 20, 21]. 

Corrections of the penguin type as shown in Figure 4 enter the perturbatively cal- 
culable hard-scattering kernels for B —)■ M1M2 decay amplitudes in QCD factorization. 
Such penguin contributions are consistently included in the factorization formula at next- 
to-leading order in perturbation theory, and at leading order in the heavy-quark limit 
[2, 3, 4]. Nonperturbative strong interactions of the intermediate quarks and gluons do 
exist, but they are suppressed by powers of Aqcd/?^?) with respect to the leading, fac- 
torizable amplitudes. In the following we shall discuss how the suppression comes about 
in the penguin contributions with a charm-quark or a light-quark loop. Of particular 
interest is the case of the charm penguin, where the validity of the usual factorization 
formula has been questioned in the literature [22]. 

In order to discuss the physics of penguin-type matrix elements with charm-loop con- 
tractions, it is instructive to consider first the case of electromagnetic penguins (Figure 4 
with the gluon replaced by a photon). These are rather similar to the more prominent 
gluonic penguins, but the QCD dynamics is simpler in this case. The electromagnetic 
penguin matrix element contributes at 0{a) to the factorization coefficient afg (p = u, 
c) in the transition operator for B irn decays, which contains the term [4] 



V-A ■ 



(45) 



The penguin contraction of the dominant current-current operators 2 gives a correc- 
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tion, included in a^Q, which reads 
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(46) 



where Sp = mi /ml refers to the quark p = u, c inside the loop, and 



with 



G-jri-s) = / G(s — ze, 1 — x) </)7r(x) 
Jo 

G(s, x) = — 4/ duu{l — u)\n\s — u{l — u)x 
Jo 

n,(o) 



(47) 



2, 

-- m s 
3 



Ug{xml) 



(48) 



is the penguin function, which involves a convolution of a perturbative hard-scattering 
kernel with the leading-twist pion distribution amplitude 4>Tr{x) [2, 3]. Eq. (46) represents 
the leading contribution of the electromagnetic penguin to the hard-scattering amplitude 
in QCD factorization, evaluated at lowest order in perturbation theory. The function 
Gt,{sp) is an 0{1) quantity, and Aafg = 0{a). Note the presence in the one-loop kernel 
G{s, x) of the vector-current correlator discussed in (9). Indeed, the expression in square 
parentheses in (46) can be rewritten as a term proportional to the convolution of the 
unsubtracted correlator Ilq{xml) with the pion distribution amplitude plus an additive 
correction of 2/3, which arises from the fact that in the NDR scheme with anticommuting 
75, which we have adopted in [4], the four- dimensional Fierz identities are violated. 

Consider now the effect of a cc resonance in the case of the charm-loop contribution. 
This corresponds to the diagram in Figure 4 (with photon exchange), but where the 
cc loop is replaced hj a ip resonance propagator. Assuming factorization at the weak 
vertex, the contribution of the intermediate resonance to the decay amplitude may be 
evaluated in a straightforward way. For the equivalent of the function G{s,x) one finds, 
using the narrow-width approximation. 



G^'i 



x) 



3m 



^ X 



+ ie 



(49) 



M^/ml. This may be integrated with the pion distribution amplitude 

(x) = Qxx one obtains (for 



where 

4>n{x) according to (47). Employing the asymptotic form 
< < 1) 



Gt 



mt 



2r,,,ri 



In- 



2r„ 



(50) 



This contribution has to be considered as part of a hadronic representation of the hard- 
scattering amplitude G-,,. It is suppressed in the heavy-quark limit. Consequently the 
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existence of the resonance cannot alter substantially the quark-level result for 6*,^, neither 
parametrically nor numerically.^ More specifically, in the Coulombic approximation the 
resonance contributes only at third order in a^, i.e. ~ a^(mc/mfe)^, which follows 
from (16) and (50). In the opposite limit, where the charm quark is treated as a light 
quark, we obtain instead ~ {-^qcd / fnt)^ ■ This demonstrates that the contribution of 
a light resonance is suppressed by two powers of A^qcd /nih. Note also that the imaginary 
part related to a light resonance has a suppression by four powers of Aqco/'^b due to 
the additional endpoint suppression from the pion distribution amplitude in (50). 

An important difference between the cases B ^ Xg and B ^ tttt results from 
the nature of the smearing procedure. In the first case the correlator is first squared 
and then integrated over phase-space, see (31). In the second case the contribution to 
the decay amplitude is given by an integral of a correlator with the pion distribution 
amplitude. This integral, for which duality holds, is then squared when one computes 
the decay rate. Note that the nature of the correlator itself is more complicated in the 
latter case. Instead of a correlation function of currents, the relevant correlator is the 
(ttttI . . . I-B) matrix element of the time-ordered product of a linear combination of the 
current-current operators Ql 2 in the effective weak Hamiltonian with Standard Model 
interactions converting the cc pair into light quarks and gluons. 

The assumption of factorization of the cc loop from the remaining quarks in the 
diagram is, of course, an idealization. However, gluon exchange between the charm loop 
and the b and s quarks does not invalidate our arguments. The most dangerous resonance 
contributions arise in the channel where the cc loop is in a colour-singlet state. Single- 
gluon exchange graphs are therefore less sensitive to such effects. In the case of the QCD 
penguin contribution the leading term in the QCD factorization approach contains the 
cc pair in a colour-octet state, and hence it is even less affected by resonance effects 
than the electroweak penguin contribution. When rricv'^ ~ Aqcd is assumed, soft gluon 
exchange between the charm loop and the b and s quarks (hadronically, the remnant of 
the decaying B meson) is expected to provide 0{1) modifications in the resonance region, 
but it does not alter the power counting. Since the contribution from the resonance region 
is small, as shown above, the factorization result at leading power in the heavy-quark 
expansion remains unaffected. 

The conclusion that the resonance contribution to the charm-penguin diagrams for 
exclusive hadronic B decays into two light mesons is parametrically suppressed in the 

''We estimate Gf « -0.18 + QMi for mt = 4.8 GeV, = 0.401 GeV. The partonic evaluation of 
the charm-loop contribution results in G^(sc) = 2.32 + 1.19i for TOc = 1.4 GeV, and in 2.29 + 1.43z 
for the smaller charm- mass value nic = 1.25 GeV. This should be compared to the result obtained 
from the hadronic representation of Il{q^) in the KS model [9] giving G^^{sc) ~ 1.86 -I- 2.02i and 
2.01 -|- 2.02z for the two charm-mass values, respectively. Thus, within uncertainties (pcrturbativc, 
parametric, power corrections), the partonic and the hadronic (KS) results agree reasonably well. Note 
that the resonance contributions are more important for the imaginary part of G7r(sc) than for its real 
part. This is consistent with the behaviour of moment sum rules to determine the charm-quark mass, 
which are dominated by the charmonium resonances even for low moment order [16]. While the above 
numbers are only estimates for the purpose of illustration, the analysis in [16] also indicates that global 
quark- hadron duality works well for Il{q^) between = and = m^. 
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heavy-quark limit, irrespective of whether the charm quark is treated as being heavy or 
hght, is in agreement with the general arguments on the suppression of the resonance or 
threshold region of cc loops given in [23] . It disagrees with [24] however, whose authors 
claim that QCD factorization fails for charm penguins even at leading power in the 
heavy-quark limit. We now explain why this claim is incorrect. The argument of [24] is 
based on the estimate 

^~a.(2m,)/(^)^, (51) 

where Acs is the non-perturbative contribution from the charm-quark penguins, v is the 
small charm-quark velocity, and / denotes some function of the ratio of heavy-quark 
masses. We find that there is an important factor missing on the right-hand side of (51). 
The matching described in the appendix of [24] applies to the threshold region and is 
therefore done for a fixed value of uml, the gluon virtuality in Figure 4. The effective 
amplitude obtained in this way must then be integrated over the threshold region of size 
~ 4m^w^, which implies an integral of the form (with r = 4m^/m^) 

duCriu)..., (52) 

r(l— v^) 

where Cf'^"'^ is the coefficient function of the cc production operator defined in [24]. 
Outside the threshold region the standard QCD factorization formulae apply and the 
factorization of charm loops is undisputed. Matching the amplitude gives Cj^°'^ ~ 1, and 
from the restricted integral in (52) we obtain a suppression factor of order 4m^f^/m^, 
by which the estimate in (51) should be multiplied to give 



a,(2mc)/^ Ux (53) 



This can easily be checked explicitly for the diagram in Figure 4. One may also note from 
(50) that the scaling of oiG^ in the non-relativistic limit is of the same form as (53). The 
reason for the omission in [24] of the last factor in (53) can be traced to the fact that the 
authors incorrectly implement the restriction to the threshold region by assuming CY'^^ — 
6{u — r) and so miss the suppression from the small size of this region. Thus, contrary 
to a statement in their paper, the "phase-space suppression factors" discussed in [23] 
are not included in (51). The analogy with inclusive quarkonium production mentioned 
in [24] is misleading, since in this case factorization applies to the cross section, which 
contains the phase-space 5-function in its definition, while here factorization is applied to 
an amplitude for fixed kinematics, which contains no (5-function. The last factor in (53) 
suppresses non-perturbative effects in the charm-penguin contribution to the amplitude. 
Following [24] and taking rricv'^ ~ Aqcd, the additional factor is readily seen to give 
a power suppression, specifically Acc/Ai^o ~ {oig/v) f{2mc/mb) {KQCj^/rribf'. If, on the 
other hand, we adopt the power counting rricv'^ ^ Aqcd, then the non-relativistic scales 
can be treated in perturbation theory, and the uncalculable non-perturbative effects are 
again power suppressed [23]. 
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We hope that this discussion clarifies that non-perturbative eff'ects in the threshold 
region in charm-penguin diagrams are power suppressed in the heavy-quark limit and 
hence do not spoil QCD factorization. 

7 Conclusions 

We have explained the origin of the well-known fact that the integrated B Xg 
branching fraction is dominated by resonance contributions from narrow charmonium 
states, such as i? ^ Xgip Xsl~^l~, which exceed the non- resonant charm- loop contri- 
bution by two orders of magnitude. Our arguments, which are based on well-established 
theoretical tools, show that quark-hadron duality cannot be expected to hold for these 
decays. On the other hand, we have shown that corresponding resonance effects lead 
to highly suppressed contributions to charm-penguin amplitudes in two-body hadronic 
B decays of the type B hit, which do not invalidate the standard picture of QCD 
factorization. 

The contributions of charmonium states to the different observables studied in this 
note are described in terms of resonance contributions to a correlation function Il{q^). 
The question of when large deviations from quark-hadron duality arise is tightly linked 
to how this correlator enters the formulae for the various decay rates. In fully inclusive 
quantities such as the cc contribution to the e~^e~ — > hadrons cross section, the imaginary 
part of a current- current correlator is integrated over phase space, and duality holds in 
the sense of [1]. A similar situation is realized in inclusive semileptonic or hadronic 
i?-meson decays, not discussed here. The inclusive decay rates are given in terms of 
forward 5-meson matrix elements of the imaginary parts of correlators of two effective 
weak Hamiltonians [25, 26]. 

In more complicated cases such as the penguin contributions to exclusive hadronic 
decays B M1M2 in the heavy-quark limit, the charm-penguin contributions to the 
decay amplitudes are given in terms of hadronic matrix elements of time-ordered products 
of current-current operators from the effective weak Hamiltonian with Standard Model 
Lagrangian insertions, convoluted with the pion distribution amplitude. Even though the 
QCD dynamics is considerably more complicated in this case, the resonance physics is 
similar, and duality still holds in a global sense, that the integrated (smeared) correlation 
functions calculated in a partonic picture approximate the true hadronic correlation 
functions in the heavy-quark limit. 

The charm-resonance contributions to B —>■ X^ decays, on the other hand, are 
expressed in terms of a phase-space integral over the absolute square of a correlator. 
We have shown that due to the fact that the widths of the lowest-lying charmonium 
states are very narrow, this integral violates the duality hypothesis by an amount that 
becomes arbitrarily large in the limit 0. Our main conclusions are derived with 

the example of a toy model, which keeps the essential features of the resonance physics 
while simplifying other aspects of the decay processes. 

Our discussion in this note has concentrated on the inclusive leptonic process B — >• 
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Xs for which duahty violations in the charm resonance region are particularly pro- 
nounced and experimentally well studied. An analogous discussion applies to the related 
radiative process B Xg'y. 
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